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High-resolution optical spectroscopy of DY Cen: diffuse 
interstellar bands in a proto-fullerene circumstellar environment? 
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ABSTRACT 

We search high- resolution and high-quality VLT/UVES optical spectra of the 
hot R Coronae Borealis (RCB) star DY Cen for electronic transitions of the Ceo 
molecule and diffuse interstellar bands (DIBs). We report the non-detection of 
the strongest Ceo electronic transitions (e.g., those at ~3760, 3980, and 4024 A). 
Absence of Ceo absorption bands may support recent laboratory results, which 
show that the ~7.0, 8.5, 17.4, and 18.8 /im emission features seen in DY Cen - and 
other similar objects with PAH-like dominated IR spectra - are attributable to 
proto-fullerenes or fullerene precursors rather than to Ceo- DIBs towards DY Cen 
are normal for its reddening; the only exception is the DIB at 6284 A (possibly 
also the 7223A DIB) that is found to be unusually strong. We also report the 
detection of a new broad (FWHM~2 A) and unidentified feature centered at 
~4000 A. We suggest that this new band may be related to the circumstellar 
proto-fullerenes seen at infrared wavelengths. 

Subject headings: astrochemistry — circumstellar matter — ISM: molecules - 
dust, extinction — Stars: chemically peculiar — Stars: individual: (DY Cen) 

1. Introduction 

Fullerenes and fullerene-related molecules have attracted much attention since their 
discovery in the laboratory (Kroto et al. 1985) due to their potential applications in su- 
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per conducting materials, optical devices, and medicine. In the astrophysical context, these 
complex molecules have been proposed as explanations for unidentified astronomical features 
such as the intense UV absorption band at 217 nm (e.g., Cataldo & Iglesias-Groth 2009) 
and the enigmatic diffuse interstellar bands (DIBs) (Herbig 1995, 2000; Iglesias-Groth 2007). 
The remarkable stability of fullerenes against intense radiation (e.g., Cataldo, Strazzulla & 
Iglesias-Groth 2009) has encouraged the idea that fullerenes should be present in the in- 
terstellar medium. Indeed, the 9577 and 9632 A DIBs observed in a few reddened stars 
lie near two electronic transitions of the C 60 cation observed in rare gas matrices (Foing 
& Ehrenfreund 1994) but the presence of fullerenes in astrophysical environments has been 
uncertain until recently, when the mid-IR spectral signatures of the C 60 and C 70 fullerenes 
were detected in several Planetary Nebulae (Cami et al. 2010; Garcia-Hernandez et al. 2010, 
2011c). 

Interestingly, Garcia-Hernandez et al. (2010) showed that contrary to theoretical and 
experimental expectation, fullerenes are efficiently formed in normal H-rich circumstellar 
environments. The prevailing laboratory view was that the synthesis of fullerenes is much 
more efficient under hydrogen-poor conditions (Kroto et al. 1985; De Vries et al. 1993). 
Thus, circumstellar envelopes around R Coronae Borealis (RCB) stars (see e.g., Lambert & 
Rao 1994 for a review) have been considered as promising environments for the formation 
of fullerene molecules (Goeres & Sedlmayr 1992). This is because the RCBs' hydrogen 
deficiency together with the He- and C-rich character of the gas resembles the laboratory 
conditions where fullerenes are produced. However, this suggestion is not supported by 
our recent Spitzer observations of a large sample of RCB stars, which show that fullerene 
formation is inefficient in the highly H-deficient environments characteristic of most RCBs 
(Garcia-Hernandez et al. 2011a,b). 

Surprisingly, the only exception among RCBs is the least H-deficient RCB star DY 
Cen (possibly also the second least H-deficient RCB star V854 Cen) that shows mid-IR 
emission features at ~7.0, 8.5, 17.4, and 18.8 /j,m coincident with the IR transitions of Ceo- 
Contrary to fullerene-containing PNe, the latter Ceo mid-IR features in DY Cen are seen in 
conjunction with very strong polycyclic aromatic hydrocarbon (PAH) features. This striking 
difference between the IR spectra of fullerene-containing PNe and DY Cen is explained by 
Duley & Hu (2012) who attribute the ~7.0, 8.5, 17.4, and 18.8 pim features in DY Cen to 
proto-fullerenes (or fullerene precursors) rather than to the Ceo molecule. Based on their 
laboratory spectroscopy of hydrogenated amorphous carbon (HAC) nanoparticles, Duley & 
Hu (2012) suggest that these four IR features may be due to proto-fullerenes in sources 
such as reflection nebulae, RCB stars, and proto-PNe, which also show the unidentified 16.4 
/iin feature and other very strong PAH-like features. In light of the now suspect infrared 
identification of neutral Ceo, we sought to detect the molecule's electronic transitions in 
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absorption from DY Cen's optical spectrum, as an additional check to the presence of C 60 
and see whether Duley & Hu (2012) 's "proto-fullerenes" suggestion could be substantiated. 

2. VLT/UVES optical spectroscopy 

High-resolution (R~30,000) and high-quality (signal-to-noise ratio S/N>200) optical 
(~3300-9450 A) spectroscopic observations of DY Cen (V=12.7, E(B-V)=0.50, see Table 1) 
were obtained in the period February-March 2010 at the European Southern Observatory 
(ESO). The UVES spectrograph at the ESO VLT (Paranal, Chile) was used with the 1.2" 
slit and the standard setting DIC2 (390+760). We obtained 11 individual exposures of 1800 
s each, giving a total exposure time of 5.5 hours. The S/N in the continuum in the summed 
DY Cen's spectrum is ~200 at 4000 A and higher than 250 at wavelengths longer than 6000 
A. The observed spectra - processed by the UVES data reduction pipeline (Ballester et al. 
2000) - were corrected for heliocentric motion and the stellar continuum was fitted by using 
standard astronomical tasks in IRAF. 

As a comparison star for DY Cen in our search for Ceo and DIBs, we selected the very 
nearby B0. 5 la bright supergiant HD 115842, which only differs in Galactic longitude and 
latitude by 0°.88 and 1°.46, respectively. This is because HD 115842 has the same reddening 
(E(B-V)=0.5) as DY Cen and it may be located at a comparable distance, sampling similar 
ISM conditions. A comparison between HD 115842 and DY Cen will probably enlighten the 
environmental changes associated with DY Cen. Thus, HD 115842 was also observed on the 
same dates as DY Cen with the same VLT/UVES set-up. Table 1 compares DY Cen with 
the comparison star HD 115842 in terms of Galactic coordinates, spectral type, magnitude, 
and reddening. Finally, a spectrum of BD —9° 4395, an extreme helium (EHe) star that 
shows a photospheric spectrum similar to that of DY Cen, was obtained by us with the 2.7m 
Harlan Smith telescope and the Tull coude spectrograph (Tull et al. 1995) at a spectral 
resolving power of 60,000. 

3. A search for neutral C 60 

We have inspected the optical spectrum of DY Cen in order to search for the presence 
of optical absorptions of neutral Ceo- The allowed transitions at wavelengths below 4100 
A are much more intense than those of the Herzberg- Teller induced forbidden transitions 
at longer wavelengths (Leach 1992). The strongest allowed electronic transitions of neutral 
gas phase C 6 o molecules, as predicted from laboratory experiments, are located at 3760 ± 5, 
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3980±0.5, and 4024.0±0.5 A with widths of 8, 6, and 4 A, respectively (Sassara et al. 2001). 
The oscillator strength (f) of the ~4024 Cqq band is similar to that at 3980 A (Leach 1992), 
while the 3760 A band has a five times higher oscillator strength (Braga et al. 1991). DY 
Cen's optical spectrum is dominated by strong He I, C II, Ne I, O II, and N II absorption 
and variable emission lines together with nebular forbidden emission lines from [O II] , [N II] , 
[S II], etc. The presence of these absorption and emission lines (FWHM~1 A) complicates 
the detection of broad and weak absorption features. 

We can find no evidence for the presence of Ceo in absorption (or emission) at the 
wavelengths of the expected electronic transitions: 

The 3760 ±5 A band: This, the strongest of the three in our bandpass, falls amongst 
a series of O II lines. These lines give consistent O abundances suggesting that none is 
does blended with a Ceo feature (Figure 1, left panel), even though its central wavelength is 
uncertain by ±5 A (Sassara et al. 2001). Note that the apparent additional absorption at 
3756 A in DY Cen is a stellar blend (Figure 1, left panel). 

The 3980±0.5 A band: This coincides with a stellar C II line from a multiplet represented 
by other lines. A secondary contributor may be the Al III 5d 2 D - 8f 2 F° multiplet (Figure 
1, right panel). The C II lines are present also in BD +9°4395. This coincidence and the 
line's 'stellar' width show that Ceo is n °t a significant contributor. An interesting feature of 
the right panel of Figure 1 is the appearance of an unidentified broad feature in DY Cen at 
4000 A (FWHM ~ 2 A; equivalent width of 91 mA). It is to be noted here that the 4000 A 
feature is real. First, a strong argument in favour of the 4000 A feature being real is that it 
doest not appear in the spectra of the comparison star taken with the same set-up and at 
the same time. Second, the 4000 A band is clearly seen when reducing the corresponding 
Echelle order alone. This is not a known DIB; it's not present in HD 115842 where all DIBs 
have a similar strength to those in DY Cen (see below) and neither is it listed by Hobbs et 
al. (2008). Identification as Ceo would demand a wavelength error of either 20 A for the 
3980 A transition or 24 A for the 4024 A transition but Sassara et al. give the wavelength 
uncertainty as just ±0.5 A for both transitions. 

The C 6 o band at 4024. 0±0. 5: This is irretrievably blended with the He I 4026 A line and 
adjacent lines. 

We estimate that the one-sigma detection limits on the column density from the 3760 
and 3980 A C 60 bands in our DY Cen spectra (S/N~200) are 4.6 x 10 12 and 1.5 x 10 13 
cm -2 , respectiveljIE It is of interest to compare this column density limit of about 10 13 cm -2 



1 l-sigma detection limits for the EQWs in our spectra scale as ~ 1.064 x FWHM / (S/N) (see e.g., Hobbs 
et al. 2008). 
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with the number of Ceo molecules calculated from the infrared emission features attributed 
to these molecules. This exercise is attempted next. 

The IR features measure the number of molecules in excited vibrational states but 
the electronic transitions sought here depend on the molecules in their ground vibrational 
state. Analysis of the strengths of the IR features with the estimated A-values (Cami et 
al. 2010; Garcia-Hernandez 2011a) shows that the excitation temperature is about 600 K. 
With this temperature, we estimate the number of molecules in the ground state: iV(C60) = 
1.0 x 10 46 d 2 where d, the distance to DY Cen, is given in kpc. 

The column density pertinent to our non-detection of Ceo depends on how the line of 
sight to the star intersects the Ceo-containing part of the circumstellar shell. Assume a 
spherical shell centered on the star with inner and outer radii of R out and Ri n , respectively, 
with the line of sight running directly from R out to R in . Assume a uniform density of Ceo 
molecules throughout the shell, then the density of molecules is n(Ceo) where 



n(C m ) = pa X _ 10 ff3 N ^o) (1) 

where R out and R in are given in AU. 

Shell radii are estimated from the run of dust temperature with distance from the star. 
The equilibrium temperature of a gray dust grain T d (R) in an optically thin circumstellar 
environment is given by 



r, \ 0.5 



= t, ( 2) 

where R* is the stellar radius, R is the radial distance from the stellar center, and T* is the 
stellar blackbody temperature (Kwok 2007, p. 314, Equation (10.32)). If stellar radii and 
temperature are estimated from the stellar luminosity L/L Q , the distance R is provided by 

for R in AU. 

Adopting L/L & = 10 4 , we find T d = 600 K at R= 21 AU and T d = 270 K at 104 AU 
where the former temperature is the excitation temperature of the Ceo molecules and the 
latter temperature is the blackbody temperature of the dust. If we adopt R in ~ 10 AU and 
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Rout ~ 100 AU, we obtain n(C 60 ) = 0.7d 2 cm -3 or n(C 60 ) = 34 cm -3 for d = 7 kpc and 
17 cm -3 for d — 5 kpc. Then, the estimates for the Ceo column density along the path 
(R out — Rin) are 4.6 x 10 16 cm~ 2 for 7 kpc and 2.3 x 10 16 cm -2 for 5 kpc. 

Such estimates are uncomfortably greater than the observed upper limit of around 
10 13 cm -2 . One possible explanation of this 1000- fold discrepancy is that the line of sight 
to DY Cen does not intersect the Cgo-containing regions of the circumstellar envelope. An 
alternative view raises doubts about the IR detection of the C 60 molecule in DY Cen (possibly 
also in V854 Cen; Garcia-Hernandez et al. 2011a), supporting the recent claim that the ~7.0, 
8.5, 17.4, and 18.8 /zm features seen in astronomical sources with PAH-like dominated IR 
spectra should be attributed to proto-fullerenes or fullerene precursors rather than to Gqq 
(Duley & Hu 2012). The latter seems to be consistent with the characteristics of the DIBs 
towards DY Cen, which are studied below. 



4. Diffuse bands from fullerene precursors? 

As we have mentioned above, we detect in DY Cen a broad unidentified feature centered 
at ~4000 A, which is seen in DY Cen only (Fig. 1). Note that no diffuse interstellar bands 
(DIBs) are known at this wavelength (see e.g., Hobbs et al. 2008). To our knowledge, no 
molecule is known to exhibit a strong electronic transition at ~4000 A and this is corrobo- 
rated by inspecting the NIST Chemistry WebBook Databas^E We speculate that this 4000 
A band may be related to the proto-fullerenes or fullerene precursors seen in the IR and that 
could be very abundant in the circumstellar envelope of DY Cen. Thus, the 4000 A feature 
may represent a new and very unusual diffuse band (DIB). 

We used the exhaustive list of DIBs provided by Hobbs et al. (2008) in HD 204827 to 
search for them in the spectrum of DY Cen. DIBs detected in DY Cen (and the comparison 
star HD 115842) are given in Table 2, where we give the central wavelength (after applying 
the heliocentric correction), full width at half- maximum (FWHM) in A, equivalent width 
(EQW) in niA, central depth (A c ), and the S/N ratio in the neighbouring continuum. Figure 
2 compares a selection of DIBs in DY Cen with those in HD 115842. Table 2 (also Figure 
2) shows that the DIBs towards DY Cen have about the same or slightly less (up to 30 per 
cent) strength, the same FWHM, and a similar radial velocity as those towards HD 115842. 
Since both stars have roughly the same interstellar reddening E(B-V), the similar profiles 
and strengths for their DIBs is not surprising. In addition, for the well studied DIBs at 
5780, 5797, 5850, 6196, 6379, and 6614 A, the normalised equivalent widths EQW/E(B-V) 



http : / / webbook. nist .gov / chemistry / 
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in DY Cen agree very well with those measured in field stars by Luna et al. (2008) and in 
the prototype star HD 183143 (Herbig 1995). 

The exceptions to the above trends are the DIBs at 6284 A and at 7223 A. The 7223A 
DIB is stronger in DY Cen than in HD 115824 as seen from the divided spectrum, but 
we could not estimate its total absorption because of the interference from the superposed 
telluric lines. Figure 3 (left panel) shows the region of 6284 A for both stars. The spectra 
have been cleaned of telluric lines. It is clear that this DIB towards DY Cen is stronger 
than towards HD 115842, suggesting that the carrier of the 6284A DIB (along with 7223 
A) is different from the rest of the classical DIBs. Indeed, the strength of this DIB seems 
to be poorly correlated with the interstellar reddening. For example, Luna et al. (2008) 
and Herbig (1995) obtained EQW/E(B-V) of 0.90 and 1.50 for field reddened stars and HD 
183143, respectively, while we measure a value of 1.10 towards DY Cen; Hobbs et al. (2008) 
give a value of 0.41 in HD 204827. Curiously, the post-AGB star IRAS 06530-0213, which 
shows the strongest unidentified 21 /im feature (Zhang, Kwok & Hrivnak 2010), also shows an 
unusually strong 6284A DIB (see Luna et al. 2008). The 21 /im feature may be attributed to 
HACs or their decomposition products such as fullerene precursors or intermediate products 
(see e.g., Garcia-Hernandez 2012). It is to be noted here that examination of atomic and 
molecular lines towards the sight lines of DY Cen and HD 115842 show that ionized lines 
are enhanced in strength in the sight line towards DY Cen relative to that of HD 115824, 
which might suggest that the carrier(s) of the 6284A and 7223A DIBs could even have some 
contribution from ionized species. 

Finally, we would like to remark that the very broad 4430 A DIB in DY Cen is almost 
identical to that in the comparison star HD 115842. This is shown in the righthand panel 
of Figure 3 where we compare the 4400—4450 A spectral region for both stars. The 4430 A 
DIB has been linked to fullerenes bigger than C 60 (e.g., C 8 o, C 2 4o, C 32 o, and C 54 o) and/or 
buckyonions such as C 6 o@C 2 4o and C 60 @C 2 4o@C 5 4o (Iglesias-Groth 2007). Our finding also 
would be consistent with fullerenes and fullerenes-containing molecules not being especially 
overabundant towards DY Cen. 



5. Concluding Remarks 

Duley & Hu (2012) show that the infrared features at ~7.0, 8.5, 17.4, and 18.8 /iin 
detected in astronomical sources with PAH-like dominated IR spectra - such as those of the 
RCB stars DY Cen and V854 Cen (Garcia-Hernandez et al. 2011a), the Reflection Nebulae 
NGC 7023 and NGC 2023 (Sellgren et al. 2010) or the proto-PN IRAS 01005+7910 (Zhang 
& Kwok 2011) should be attributed to fullerene precursors (or proto-fullerenes) rather than 
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to Ceo- Our non detection of neutral Ceo in DY Cen may support this claim . Thus, the 
carrier of the new 4000 A absorption band detected in DY Cen may be intimately related 
with the fullerene precursors seen at infrared wavelengths. Duley & Hu (2012) suggest 
that fullerene precursors contain pentagonal carbon rings; the infrared emission features 
at ~7.0 and 16.4 /iin usually detected in objects with PAH-like dominated IR spectra are 
characteristics of pentagonal rings (Moutou et al. 2000). Pentagonal carbon rings are present 
in carbon nanoparticles (e.g., HACs; Duley & Hu 2012) and nanotubes although these have 
less perfect structures than C 60 and other fullerenes. 

In summary, the non detection of neutral Ceo i n the high-quality VLT/UVES DY Cen's 
spectrum may support recent experimental work, showing that the ~7.0, 8.5, 17.4, and 
18.8 jum IR features seen in sources with PAH-like dominated spectra have to attributed 
to proto-fullerenes rather than to neutral C 60 . In addition, the new 4000 A DIB reported 
here (possibly also the carriers of the classical 6284A and 7223A DIBs) may be related to 
fullerene precursors; an organic compound containing pentagonal rings. These pentagonal 
carbon rings are usually present in HAC nanoparticles and nanotubes, suggesting that they 
may be intimately related with the formation process of fullerenes. 

This work is based on observations obtained with the ESO Very Large Telescope at 
Paranal under the Director's Discretionary Time (DDT) programme 284. D- 5048(A) (PI: D. 
A. Garcia-Hernandez). D.A.G.H. acknowledges support provided by the Spanish Ministry 
of Economy and Competitiveness under grant AYA— 2011— 27754. D.A.G.H. also thanks his 
son Mateo for his great patience during the realization of this work. D.L.L. also wishes to 
thank the Robert A. Welch Foundation of Houston, Texas for support through grant F-634. 

Facilities: VLT:UVES. 
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Table i: Observational parameters of DY Cen. 



Star 


1 


b 


V 


B-V 


E(B-V) 


Sp. Type 


My 


Distance 
(pc) 


Ref 


DY Cen 


307.958 


8.293 


12.7 


0.33 


0.47-0.50 


BO la 




4800 


1 


HD115842 


307.080 


6.834 


6.04 


0.30 


0.50 


B0.5 la 


-8.0 


3581 


















-6.9 


2157 


2 



1:- De Marco et al(2002), 2: Krelowski et al (2010), Hunter 
et al (2006) 
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Table 2: Diffuse Interstellar Bands in DY Cen and HP 115842" 



DY Cen 










HD 115842 












FWHM EQW 


Ac 


S/N 


Ac 


FWHM 


EQW 


Ac 


S/N 


(A) 


(A) 


(mA) 






(A) 


(A) 


(mA) 






4000.38 


1.92 


91 


0.04 


188 










412 


4429.66 


21.10 


1034 


0.05 


273 


4428.50 


23.00 


932 


0.04 


436 


5780.46 


2.14 


256 


0.10 


209 


5780.39 


2.10 


257 


0.11 


307 


5796.92 


0.91 


72 


0.07 


199 


5796.90 


0.89 


105 


0.11 


264 


5849.67 


0.82 


32 


0.02 


272 


5849.58 


0.83 


43 


0.05 


330 


6065.11 


0.61 


6 


0.01 


312 


6065.00 


0.64 


7 


0.01 


635 


6089.62 


0.53 


7 


0.01 


282 


6089.59 


0.53 


13 


0.02 


490 


6112.7 


1.00 


12 


0.01 


300 


6112.97 


0.62 


12 


0.02 


490 












6116.62 


0.94 


8 


0.01 


490 


6139.75 


0.82 


5 


0.01 


267 


6139.72 


0.59 


8 


0.01 


570 


6195.80 


0.60 


27 


0.04 


284 


6195.80 


0.50 


33 


0.06 


636 


6202.81 


1.28 


62 


0.04 


336 


6202.83 


1.35 


64 


0.04 


636 


6233.79 


0.79 


10 


0.01 


336 


6233.80 


0.66 


14 


0.02 


560 


6269.70 


1.21 


36 


0.02 


276 


6269.80 


1.38 


70 


0.04 


582 


6283.53 


4.29 


543 


0.12 


276 


6283.75 


4.43 


413 


0.09 


582 


6375.97 


0.92 


13 


0.02 


273 


6375.91 


0.90 


29 


0.03 


340 


6379.25 fe 


1.04 


53 






6379.04 


0.64 


70 


0.10 


340 


6597.26 


0.53 


8 


0.01 


295 


6597.19 


0.58 


9 


0.01 


378 


6613.41 


1.03 


83 


0.08 


295 


6613.45 


1.04 


149 


0.14 


304 


6660.60 


0.64 


12 


0.02 


314 


6660.53 


0.80 


39 


0.05 


265 


7223.75 




















8620.88 


3.75 


95 


0.02 


106 


8620.94 


3.92 


105 


0.02 


243 


a The 3-ct 


erros in 


the EQWs 


scale like ~3 


x FWHM / (S/N) while we 







estimate that the FHWMs are precise to the 0.03 A level; b blended 
with N II 6379.98 A 



13 




3740 



3750 3760 
X (A) 



3770 



3780 



DYCen He I 

HD1 15842 

BD-9 4395 



3990 4000 4010 

A (A) 



Fig. 1. — The spectra of DY Cen (in red) and HD 115842 (in blue) around 3760 A (left 
panel) and 4000 A (right panel; where the eHe star BD —9° 4395 is also displayed in green) 
are shown. The expected positions of the Cgo features are marked on top of the spectra.. 
The FWHM of the features are indicated by the dotted line. Note that there is no evidence 
(additional absorption) in DY Cen for the presence of the neutral Cgo features at 3760, 3980, 
and 4024 A. However, there is an additional absorption band at 4000 A in DY Cen (marked 
by a dark line and written DIB4000), which is not present in either HD 115842 or in the 
extreme He star BD —9° 4395 that is expected to show similar spectrum as DY Cen. 
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Fig. 2. — Profiles of six DIBs of various strengths and widths are plotted with respect to 
the LSR radial velocity. Note that the profiles are shifted vertically for clarity. The [Fe II] 
profile at the bottom shows the DY Cen's systemic radial velocity. Note that the strengths 
and shapes of 5780 and 6203 A DIBs are almost identical for both stars. The other profiles 
show that DIBs towards the DY Cen sight of line are weaker than towards HD 115842. To 
convert the velocity scale to heliocentric, 4.6 and 5.1 km s _1 should be added to the LSR 
velocities of DY Cen and HD 115842, respectively. 
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Fig. 3. — Spectral regions around the DIBs at 6284 A (left panel) and 4428 A (right panel) 
in DY Cen (in red) and HD 115842 (in blue). The telluric line corrected spectrum of both 
stars are shown by smooth lines (red and blue). The not fully telluric corrected spectrum 
with emission of [O I] line in DY Cen is shown in green. It is obvious that the 6284 A DIB 
is stronger in DY Cen spectrum than towards HD 115842 while the DIB at 4428 A (right 
panel) is of similar strength in both stars. 



